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A SIMPLE  APPROXIMATION  FOR  THE  FUNDAMENTAL  FREQUENCIES  OF 
TWO-SPAN  AND  THREE-SPAN  CONTINUOUS  BEAMS 

by 

A.  S-  Veletsos  and  N.  M.  Newmark 
SYNOPSIS 

A rapid  approximate  method  i6  presented  for  calculating  the 
fundamental  frequencies  of  flexural  vibration  of  two-span  beams  and  of 
particular  arrangements  of  three -span  beams  which  are  continuous  over 
non -deflecting  supports  and  are  elastically  restrained  against  rotation 
at  their  end  supports.  The  end  restraints  may  be  provided  by  actual 
coll  springs  or  they  may  represent  the  effect  of  adjoining  members , but 
in  all  cases  the  stiffnesses  of  these  restraints  are  assumed  to  be  posi- 
tive. The  mass  per  unit  of  length  and  the  flexural  rigidity  of  the  beams 
may  vary  from  one  span  to  the  next,  but  in  any  one  span  these  quantities 
are  considered  constant.  Two  numerical  examples  are  included  to  illustrate 
the  application  of  the  method. 

SIGN  CONVENTION 

The  following  sign  convention  is  used.  Clockwise  rotations  are 
taken  as  positive.  Bending  moments  at  the  ends  of  a span  are  considered 
positive  when  acting  in  a clockwise  direction  on  the  beam. 


2. 


BASIS  AMD  GENERAL  DESCRIPTION  OF  METHOD 


When  a continuous  beam  is  in  a state  of  free  oscillations,  each 
of  the  spans  is  elastically  restrained  against  rotation  at  its  ends  by 
the  rigidity  of  the  contiguous  spans  and  vibrates  with  the  same  frequency 
as  that  of  the  continuous  system.  Therefore,  the  problem  of  determining 
the  natural  frequencies  of  a continuous  beam  is  basically  the  same  as  that 
of  determining  the  corresponding  frequencies  of  one  of  its  spans  only, 
with  proper  consideration  of  the  actual  restraints  existing  at  its  ends. 
The  stiffnesses  of  these  restraints  depend  on  the  properties  of  all  the 
spans  and  on  the  order  of  the  desired  natural  frequency. 

Consider  a continuous  beam  oscillating  in  its  fundamental  mode 
of  free  vibration.  Let  the  supports  be  numbered  consecutively  starting 
with  1 at  one  end  and  terminating  with  z at  the  other  end.  Let  0^  be 
the  rotation  of  the  beam  at  an  interior  support  J,  and  M.  . , and  M,  . 
be  the  internal  bending  moments  at  end  ^ of  the  span  between  (j-l)  and 
J,  and  that  between  J and  (j+l) , respectively.  The  relationship  between 
these  quantities  may  be  expressed  by  the  equations 

- -Ki,J-l63  “d  * -KJ,j+lV  (1) 

in  which  K.  , . and  K.  are  the  stiffnesses  of  the  internal  restraints 

j,j-x  OjJ+x 

at  support  j_.  For  a hinged  condition  K = 0,  whereas  for  a fixed  con- 
dition K = infinity.  The  negative  signs  in  the  foregoing  expressions 
denote  that  for  a positive  restraint  (positive  value  of  K),  the  moment 
exerted  by  the  restraint  on  the  span  acts  in  a direction  opposite  to 


the  direction  of  rotation  of  the  span.  The  end  moments  M and  M 

1,2  z,z-i 

are  related  to  the  end  rotations  9 and  6 by  expressions  similar  to 

x z 


3. 


those  given  in  Eq.  (l).  It  will  be  assumed  that  the  stiffnesses  of  the 

end  restraints,  K and  K ,>  are  positive  and  known. 

1,2  z,z-i 

For  a natural  mode  of  free  vibration,  no  external  moment  acts 
on  the  system;  therefore, 


+ 1 " 0 

(2) 

i,d-i  + KJ, J+i  “ 0 * 

(3) 

whence 


Expressed  in  words,  Eq.  (5)  states  that  the  sum  of  the  stiffnesses  at  a 
Joint  is  equal  to  zero.  It  should  be  pointed  out  that  this  relationship 
holds  true  not  only  for  the  fundamental  mode,  but  for  the  higher  natural 
modes  as  well. 

The  procedure  to  be  presented  consists  of:  (a)  isolating  from 

the  continuous  beam  the  span  from  J to  (j+l)  subjected  to  positive  end 

restraints:  (b)  determining  the  stiffnesses  of  these  restraints,  Kj  j+1 

and  K. and  (c)  evaluating  the  fundamental  frequency  of  the  continuous 

j+x,j 

beam  from  the  approximation'1' 


f « 


r i , i i r i , 1 i 

L 2 5 + 0.  JLA  + 2 5+  B,.,  , J 


1 


J+1,J 


2Lj 


\ 


E.I 


(V 


in  which  L,,  E.I.,  and  m.  are,  respectively,  the  length,  the  flexural 
J J J J 

rigidity  of  the  cross  section,  and  the  mass  per  unit  of  length  of  the 
span  between  J and  (j+l),  and  0,  . . and  0 , are  dimensionless 

quantities  related  to  the  stiffnesses  of  the  end  restraints  by  the 
equations 


1.  "A  Simple  Approximation  for  the  Natural  Frequencies  of  Partly 
Restrained  Bars,"  by  R.  M.  Newmark  and  A.  S.  Veletsos,  Journal 
of  Applied  Mechanics,  Vol.  19,  1952,  p.  563- 
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J,j+1 


K 


J,J+!  Ejlj 


and 


0 


J+lrJ 


>1,J  Ejlj 


(5) 


The  frequency  f is  expressed  in  cycles  per  second.  Eq.  (4)  is  applicable 
to  positive  restraints  only;  it  is  for  this  reason  that  the  isolated  span 
must  be  positively  restrained. 


TWO-SPAN  BEAMS 


For  a two -span  beam,  such  as  that  shown  in  Fig.  1,  it  is  only 

necessary  to  determine  the  stiffness  of  the  restraint  exerted  by  one 

span  upon  the  other.  Let  f^  and  f^  be  the  fundamental  frequencies  of 

spans  (1,2)  and  (2,3),  assuming  that  the  beam  is  hinged  at  support  2 

({3-  ..  ip,’  0).  These  frequencies  may  readily  be  evaluated  from 
2,1  2,3 

Eq.  (4). 


If  the  supports  are  numbered  so  that  f 2 — f^,  the  stiffness 

K of  the  restraint  exerted  by  the  dynamically  stiffer  span  (l,2)  on 
2,3 

the  dynamically  weaker  span  (2,3)  will  be  greater  than  or  equal  to  zero, 
and  the  fundamental  frequency  f of  the  continuous  beam  will  lie  between 

f2  011(1  fl 

From  the  results  of  numerical  calculations  based  on  exact 
solutions,  the  following  empirical  approximation  has  been  found  for 


2,2 3 


= (K. 


t 1 ' - 


(6) 


in  which  (K0  ,)  is  the  stiffness  of  the  restraint  provided  by  span 

dyj  & 

(l,2)  under  static  conditions.  It  can  readily  be  shown  that 


2.  See  for  example  "A  Direct  Method  of  Moment  Distribution,"  by  T.  Y. 

Lin,  Transactions  A.S.C.E.,  Vol.  102,  1937,  P*  565. 
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With  K,  known  and  K , determined  from  Eq.  (6),  span  (2,3) 

2,3 

may  now  he  treated  as  a har  subjected  to  positive  end  restraints,  and 
its  fundamental  frequency,  which  is  also  the  desired  frequency  of  the 
continuous  beam,  may  be  evaluated  from  Eq.  (4)  In  this  case  j = 2 and 
J+l  = z = 3. 

The  accuracy  of  Eq,  (6)  and  that  of  the  natural  frequencies 
determined  by  the  foregoing  procedure  have  been  checked  for  over  three 
hundred  representative  beams  having  end  restraints  in  the  range  between 
hinged  and  fixed  conditions  and  spans  with  ratios  of  lengths , ratios  of 
flexural  rigidities  of  cross  section,  and  ratios  of  masses  per  unit  of 
length  in  the  range  between  zero  and  one.  The  greatest  error  was  found 
to  occur  in  the  case  of  beams  which  have  a ratio  of  flexural  rigidities 
of  cross  section  from  about  0.2  to  0.4  and  have  the  extreme  end  of  the 
dynamically  stiffer  span  hinged  or  practically  unrestrained  and  the  end 
of  the  other  span  clamped  or  very  nearly  fixed 

As  an  indication  of  the  accuracy  of  Eq.  (6)  some  representative 
results,  including  those  for  which  the  error  is  maximum,  are  given  in 
Fig.  2 In  this  figure,  the  abscissas  K / (K  ) were  determined 

2,3  2,3  s 

from  the  exact  solution,  whereas  the  quantities  f^  and  fg  for  the  ordin- 
ates were  computed  from  Eq.  (4).  The  vertical  distances  between  the 
various  points  in  this  figure  and  the  diagonal  line  represent  the  error 
involved  in  Eq.  (6).  These  particular  results  are  applicable  to  two- 
span  beams  simply  supported  at  one  end  and  elastically  restrained  at  the 
other.  It  should  be  noted  that  for  the  limiting  values  of  fg/f^  = 0 and 


6. 


fg/f^  * 1.00  Eq.  (6)  is  exact. 

Figure  2 indicates  that,  when  the  ratio  of  the  flexural 
rigidities  is  the  variable,  the  error  in  Eq.  (6)  is  appreciable.  How- 
ever, because  the  natural  frequencies  of  elastically  restrained  bars 
are  not  very  sensitive  to  the  stiffnesses  of  the  end  restraints,  the 
error  in  the  natural  frequencies  determined  by  using  Eq.  (6)  is  for  all 
practical  purposes  insignificant.  By  comparing  the  exact  natural  fre- 
quencies of  the  more  than  three  hundred  beams  referred  to  previously 
with  those  determined  by  the  foregoing  procedure,  it  was  found  that  the 
maximum  error  in  the  frequencies  determined  by  the  approximate  method 
is  within  ^ 5 percent . 

Example.  - As  an  illustration,  consider  a beam  having  the 
following  characteristics: 

L = 0.80L  , El  * E I , m = 0.8lm  , 

1 2 112  2 1 2 

K = 1.0E  I /L  and  K = 5-OE  I /L  . 

1,2  111  3,2  222 

The  frequencies  f and  f , determined  from  Eq.  (4) , are 
2 1 

f = 1.00  x 1.25  f = 1.25  f , 

2 0 0 

\ ^ ^ * 0^75  fo  - 1M  V 


where 


f 

o 


The  static  stiffness  of  the  restraint  exerted  by  span  (1,2) 


on  span  (2,5)  is 


(K  ) = 0.80  x 4.0  E I /L  - 5-2  E I /L  , 

2,3  S 1 l'  1 1 1 
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and  the  corresponding  dynamic  stiffness,  computed  from  Eq.  (6),  is 


Then, 


K - 0.558  x 3.2  E I /L  » 1.79  El/L. 

2,3  111  1 1 1 


P2,3  = 1,79  XO0S  2*24,  P3,2  “ 5>°» 


f = 1.155  X 1.25  f - 1.44  f . 

o o 

The  exact  value  of  f,  neglecting  the  effects  of  damping, 
rotatory  inertia,  and  shearing  deformation,  is  1.43  f • 


THREE -SPAN  BEAMS 


Consider  the  three-span  beam  shown  in  Fig.  3-  Let  f°,  f°, 

1 2 

and  f be,  respectively,  the  fundamental  frequencies  of  spans  (1,2 ), 

(2,3),  and  (3,4),  assuming  that  the  beam  is  hinged  over  its  interior 

supports  (p  = P = p = P =0).  These  frequencies  are  de- 
2,1  2,3  3,2  3,4 

termined  from  Eq.  (4)  Only  those  cases  will  here  be  considered  for 
which  f°  and  f°  are  sufficiently  larger  them  f°  so  that,  when  the  beam 
vibrates  in  its  fundamental  mode,  the  restraints  exerted  on  the  central 
span  erne  positive. 

The  stiffnesses  K and  K are  determined  by  successive 

2,3  3,2 

approximations  as  follows:  One  assumes  a value  for,  say,  K and,  by 

3 ;2 

treating  the  portion  of  the  beam  between  supports  1 and  3 as  a two-span 

continuous  beam  in  the  manner  described  previously,  calculates  an 

approximate  value  for  K . Using  this  value  of  K and  working  with 

2,3  2,3 

the  portion  of  the  beam  between  supports  2 and  4,  one  then  computes  a 

new  value  for  K • From  this  revised  value  of  K , one  then  obtains 
3,2  3,2 
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a new  value  of  K • This  procedure  is  repeated  until  the  values  of  both 

2,3 

K and  K converge.  Reasonable  convergence  is  generally  obtained  in 

2,3  3,2 

two  or  three  cycles. 

Having  K and  K , the  fundamental  frequency  f of  the  con- 
2,3  3,2 

tinuous  beam  may  be  calculated  from  Eq.  (4)  by  considering  the  central 
span  as  an  elastically  restrained  bar.  As  before,  by  comparing  the  exact 
and  the  approximate  natural  frequencies  for  a number  of  representative 
beams  covering  the  possible  range  of  variables,  it  has  been  concluded 
that  the  maximum  error  in  the  value  of  f determined  by  the  foregoing  pro- 
cedure is  of  the  order  of  t 5 percent. 

Example.  - As  an  illustration,  consider  a beam  having  the 
following  characteristics. 

E I a 0.80E  I , La  O.85L  , m a 0.80m  , 
lx  221  2 1 2 


E I a 0.80E  I , La  0.9QL  , 

3 3 2 2 3 2 


m = 0.70m  , 
3 2 


K a 4.0E  I /L  and 
1,2  xxx 


= 1.6E  I /L  . 

3 3 3 


The  frequency  f ° = ~~  , E I /m  * f . 


2L“ 


222 


The  frequencies  f and  f°,  evaluated  from  Eq.  (4),  are 
x 3 

f°  = 1.22  X 1.00  X 1.584  f a 1.69  f , 
1 o x o' 


f = 1.00  X 1.12  X 1.520  f a 1.48  f . 

3 OO 


In  this  particular  case,  the  successive  approximation  pro- 
cedure is  started  by  taking  for  the  dynamic  stiffness  K a value 

3,2 

equal  to  one -half  the  corresponding  static  stiffness  (K  ) . The 

3,2  a 

value  of  the  latter  is  determined  from  Eq.  (7)  "by  replacing  the 


V'i- 


9- 


quantities  E I and  L by  E I and  L , and  0 by  0 , 

11  133  3 1,2  4,3 


(K  ) « 0.8214-  x io  E I /L  » 3-286  E I /L  - 2.92  E I /L  . 

3 j 2 S 333  333  2 £ 2 


Hence, 


K * 0.5  x 2.92  El  A * 1.46  E I /L 

1 3,2  2 2 2 2 2 


In  this  expression  K denotes  the  first  approximation  to  K .In 

13,2  3,2 

general,  K,  . . will  designate  the  value  of  K.  . . at  the  beginning  of 

the  n-th  cycle  of  the  procedure. 

The  portion  of  the  beam,  between  supports  1 and  3 is  now  treated 

as  a two-span  continuous  beam  with  K equal  to  1.46E  I /L  . The  fre- 

3,2  2 2 2 

quencies  f and  f of  the  individual  spans  (assuming  the  beam  hinged  at 
12 

support  2)  are 

f = f 0 = I.69  f , 

11  ° 

and  f * 1.11  f . 

2 0 

The  static  stiffness  of  the  restraint  provided  by  span  (1,2)  on  span 
(2,3)  is  determined  from  Eq.  (7)  as 


(K  ) - 0.875  X 4.0  E I/L 

2.3  s 111 


3-50  E I A = 3.29  E I A . 

Ill  222 


The  first  approximation  to. the  corresponding  dynamic  stiffness  is 

obtained  from  Eq.  (6)  as 

K « O.569  x 3-29  E I A ■ 1.87  El  A • 

1 2,3  222  222 


Next,  the  portion  of  the  beam  between  supports  2 and  4-  is 

considered,  with  K taken  equal  to  K . On  the  assumption  that  the 

2,3  1 2,3 

beam  is  hinged  at  support  3;  the  fundamental  frequencies  of  the  individual 
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f 


spans  are 

f = 1.136  f and  f • f°  = 1.48  f . 

2 0 3 3 o 

The  dynamic  stiffness  ^ is  obtained  from  Eq.  (6)  by  substituting 

(K  ) for  (K  ) and  f for  f , 

3,2  S 2;3  S 3 1 


K = 0.411  x 2.92  El/L  « 1.20  El/L  . 

2 3,2  222  222 

This  newly  computed  value  of  K leads  to  K = 1.91  E I /L 

3,2  2 2,3  22  2 

which,  in  turn,  leads  to  K = 1.19  E I /L  . It  should  be  observed 

3 3,2  222 

that,  for  all  practical  purposes , K is  equal  to  K and  K 

2 2,3  1 2,3  3 3,2 

is  equal  to  K . Therefore,  the  0 values  for  the  central  span  may  be 

2 3,2 

taken  as 

0 «.  1.91  and  0 * 1.19. 

2,3  3,2 

The  fundamental  frequency  f of  the  continuous  beam  is  finally 
evaluated  from  Eq.  (4),  where  j = 2,  as  follows: 

f = 1.138  x 1.096  f = 1.25  f . 

The  exact  value  of  f,  neglecting  the  effects  of  damping,  rotatory  inertia. 


and  shearing  distortion,  is  also  equal  to  1.25  f • 
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